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Abstract 

The potential of very-long-baseline interferometry (VLBI) 
is examined for use in the determination of wind speeds in 
Venus' lower atmosphere via the differential tracking of entry 
probes. A simplified mathematical model is presented in 
detail. An incomplete error analysis based on this model 
permits an educated guess to be made: An uncertainty in wind 

speed determination of no more than about 100t -1 m/sec, where 
t>l is the corresponding time resolution in seconds , is an 
achievable goal — without the use of transponders on the 
miniprobes. Certain important issues raised in the report 
must be resolved before firm conclusions can be drawn. However, 
if transponders are available on all probes, there should be 
little difficulty in estimating wind speeds with useful 
precision. 
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I. Introduction 

Can the Venus sntry probes be tracked via differ- 
ential very-long-baseline interferometry (VLBI) with 
sufficient accuracy to yield useful estimates of the 
wind speeds in Venus' lower atmosphere? A definitive 
answer cannot yet be given. The current status of the 
analysis is summarized in this report with special em- 
phasis on the main areas of uncertainty. Section II con- 
tains a brief description of the basic method, with a 
mathematical model and first-order error analysis being 
developed in Sections III and IV, respectively. The main 
conclusions and the requirements for additional analysis 
are presented in Section V, 
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II. Differential VLB I 

The VLBI technique has been used successfully for 
the past five years primarily to study the structures and 
positions of compact extragalactic radio sources. For a 
strong source/ the main limitations on the accuracy achievable 
in determining the direction to the source the result of 
(i) instabilities in the frequency standard used at the ob- 
serving sites, and (ii) phase fluctuations of the signals 
introduced by the propagation medium, mainly the earth's 
atmosphere and ionosphere. If several objects in nearly 
the same direction are observed simultaneously, these error 
sources can either be eliminated or drastically reduced in 
their effect on determinations of relative position. The 
frequency standard need only be sufficiently stable to 
allow fringes to be obtained on the strongest source which 
then acts as the standard for comparison with the signals 
from the other sources observed simultaneously. The prop- 
agation medium effects cancel to the extent that the sig- 
nals from the different sources received at a given site pass 
through identical paths in the earth's atmosphere and ion- 
osphere. Thus the accuracy in relative position determin- 
ation can exceed that of "absolute" position determination 
by several orders of magnitude. 
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In observing the quasar 3C279, for example, our VLBI 
group (Whitney et al. 1971) discovered that its structure 
was consistent with a two-point-source model; these two 
"points” were separated by about 1.5x10”^ arcseconds and 
the standard error in the determination of the separation 
was only 6x10 ^ arcseconds in the right ascension compo- 
nent. This extremely small error in relative position de- 
termination is meaningful because of the small angular sep- 
aration of the putative two point sources and the consequent 
high order of cancellation of the propagation medium effects. 

We call this technique of relative position determin- 
ation differential VLBI. We have also applied it successfully 
to the Apollo 16 Lunar Rover whose position relative to the 
Lunar Module was monitored throughout the first EVA by use 
of this method (Shapiro et ad. 1972; Counselman et al. 1972) 
Although the tracking systems were far from optimally arranged 
for the task and although the radio frequency of the Rover 
differed from that of the Module by lZtftfj, the final po- 
sition of the Rover calculated via the differential VLBI 
technique differed from the estimates of the astronauts by 
less than 30 m. No accurate intermediate check-points are 
available for comparison. 

The application of differential VLBI to the tracking 
of the Venus entry probes differs in several important 



respects from the Rover-Module case. On the positive 
side of the ledger, we have the possibilities (i) to 
choose nearly identical entry probe transmitter freq- 
uencies ( Ai £ S'O.JcH^) to insure that if the different 
signals pass through the same plasma environment, the 
latter s effect on phase path will cancel upon differ- 
encing; (ii) to design the receiver equipment so that 
at a given site the local-oscillator signals introduce 
the same phase noise when mixing with each of the probe 
signals (the commonality implies that this source of 
noise will also cancel upon differencing) ? (iii) to 
utilize a phase-coherent transponder on at least one 
of the entry probes; and (iv) to select earth-tracking 
sites with greater east-west and north-south baseline 
components. The negative side of the ledger contains 
more entries: (i) Venus will be about 200 times further 

away than the moon., causing a corresponding reduction in 
the accuracy of determination of the projected distance 
between tracked objects; (ii) the interplanetary medium 
has a much greater influence on Venus-earth than on moon- 
earth signals; (iii) Venus has an ionosphere and a thick 
atmosphere; the moon has virtually none of either; 

(iv) the tracked Venus probes will move relatively 
unconstrained through a fluid; the Rover was constrained 
to adhere to the lunar surface and hence the intrinsic 
two-dimensional differential VLBI tracking result could 
be converted to three-dimensional relative position by 
use of lunar topographic data; (v) the separation between 
Rover and Module was known at the 
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start of the VLBI tracking period? for the entry probes 
corresponding information will probably not be available 
and thus in the latter case only the monitoring of changes 
in the (projected) separations of the probes, i,e. only 
the monitoring of (projected) velocity differences, will 
be possible; andl (vi) f the thermal .environment of the entry 
probes will be far less stable than for the Rover and 
Module thus tending to cause greater variations in the 
transmitter frequencies of the probes. 

How does the differential VLBI technique compare 
with the straightforward use of a turnaround transponder? 
In fact, they are complementary: the transponder supplies 

the velocity and VLBI the transverse components of 

the velocity,* The VLBI approach can be used with either a 
transponder or a ^-running oscillator to determine the 
transverse components; the radial component cannot be use- 
fully inferred without a transponder unless the a priori 
knowledge of the transmitter frequency is sufficiently ac- 
curate. For signal propagation in a vacuum , the trans- 
ponder can have an enormous advantage:- all other aspects 
being equivalent, the error in the determination of radial 
velocity will be less than for the differential VLBI determ 


*By "radial" we mean parallel to the earth -Venus line 
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ination of the transverse components by the ratio of 
the VLBI baseline to the distance from the earth to the 
( i * e • r by the parallax) * For the Venus probes, 
this enormous advantage — approximately 2xl0 4 in accuracy- 
is offset to a great extent by the systematic errors in- 
troduced by the propagation medium which largely cancel 
in the differential VLBI procedure. One further point 
needs to be made here: the differential VLBI procedure 

which is needed to cancel these errors yields only the 
relative transverse components of velocities for a pair 
(or more) of probes; the transponder approach yields the 
"absolute" radial velocity for each probe. 


We may now address briefly the main problem — the 
determination of the wind speeds in Venus ' lower 
atmosphere. We distinguish two cases: 

W Transp onders Available on Entry Probes . Here 
we would be able to estimate usefully the 
velocity vector for all probes from the 
ordinary Doppler data. The a priori knowledge 
of both the geometry of entry and the terminal 
vertical velocity for each probe will most 
. likely be of sufficient accuracy for this 
purpose. The (two-way) effects of the 
atmosphere of Venus will introduce uncertainties 
well below the meter-per-second level 

unless the geometry is particularly unfavorable. 
The wind speeds will be given by the projection 
of the velocity vector on the plane normal to 
the local vertical at the probe's position. 

(We assume that the probe has 


‘•“ear'herl "formi r- 7 " 
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velocity in both the vertical and horizontal 
directions,) Under these circumstances, the 
VLBI measurements may not be competitive. But 
they will still be of interest to provide a 
check. 

(ii) Transponders N'ot Available in Entry Probes . 

Here essentially only the pairwise differences 
i-n the probes' velocities projected on the 
plane normal to the earth-Venus line will be 
available (except in the unlikely event that 
the transmitter frequency of one or more of the 
probes is known very accurately). In general, 
there is difficulty in separating the contrib- 
utions of the horizontal velocity components 
from the vertical components in the projections 
of the differences. If one of the tracked 
objects were following a ballistic trajectory 
(e.g., the bus on a flyby trajectory), then the 
contribution of the horizontal velocity component 
of each entry probe could be distinguished. 

Also, if one of the probes were directed towards 
the subearth point on Venus, it would be possible 
to identify part of the velocity projections 
as being due to winds. If neither of these 
: conditions applies, it appears that models of 
the terminal descent and appropriate filtering 
would be required to extract estimates of the 
wind speeds. Whether such estimates, would be 
useful has not yet been established. 
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Prom this qualitative introduction to the problems of 
the determination of the lower atmospheric wind speeds 
using differential VLB I*, we proceed in. the next section 
to the development of an appropriate mathematical model. 
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III. Mathematical Model 

Our goal here is to develop an algorithm for the de- 
termination of wind speeds in Venus 1 lower atmosphere from 
VLBI tracking data. Let us begin with the definitions of 
the relevant geometric- quantities. The vector distance 
R ip f tom the earth-tracking station i to the entry probe 
£ can be expressed in a geocentric reference frame as: 


R: p It, +p> ~r c (-e), ( 1) 

where * ev i s the vector distance from the center of mass 
of the earth to the center of mass of Venus; is the 

vector from the center of mass of Venus to the ftk probe, 
r i is the vector from the center of mass of the earth to 
the Ltk tracking station, t is the time of reception of 
the signal from the >iU probe at the cth tracking station, 
and is the time delay between the transmission of a 

signal from the |>tt probe and its reception at the iU track 
ing station. 


For the purposes of this section, we shall assume that 
the signals propagate in vacuum? in Section IV we will con- 
sider the medium effects explicitly. Thus, in the 
of vacuum propagation, the phase delay r^(c) may be found 
iteratively by means of a simple algorithm: 

■y*' = is V'w , u> 

( 3 ) 
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and 

Since the velocities involved are only of the order of 

■.4 

10 c, where c is the speed of light, one or two iter- 
ations will be sufficient to obtain the needed accuracy. 

If we assume continuous reception of signals start- 
ing from t=0, then the phase (r) of the signal received 
at station _i from probe £ may be written as 

t- tcp it) 

^ It) - f iTrfklx) cU , t-f) 

-y») 

where (-6) is the frequency transmitted by probe p at time 
t. From the measurements <}>. ({) we wish to estimate- the wind 
speeds, but in such a manner that we cancel to as high a de- 
gree as possible the adverse effects of the propagation med- 
ium (which is, however, ignored in the explicit formulation 
given in this section) . By the formation of symmetric double 
differences, we can insure the tendency to cancel of any po- 
tential source of error that is common either to all receiv- 
ers or to all transmitters. 

Before applying this principle, we must consider a 
means for improvement of the estimate of (t) which is not 
known accurately a priori . Our results will turn out to 
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be relatively insensitive to those estimates and so we 
will employ sums of the different (c ) 

to determine each Since f^ will vary with time in 

an unknown manner, we will estimate an average value 
^ ^ \?y appropriate for each time resolution interval of 
interest. Without any important loss in generality, we 
may consider this interval to be a constant, T (see 
Stct'io/I IE’..) . Thus, we may use Eq. (5) to obtain, 
successively : 

<fc P (t„) - (£*-.') 


h 


-^(W) V^O , (<) 


where 


= »T ; w = 

0, i, “ • 

* 

(7) 

z + T = " 

r 

a 

> 

(?) 


lo-wM) 




c: i r if -f , C/oj 
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and where ^ signifies the time derivative of the 

phase delay evaluated at t=t n _ 1/2 . if the total number 

tracking stations is I, then we consider 


K ( t, ( h) > 


1 | (t h „, tJ 

* it ' 



to be the average value of the transmitter frequency f 
over the corresponding receiving time interval 
on the assumption that the error in the measurement 
is independent of dL. If there is a dependence, a more 
suitable weighting function can easily be substituted. 


Since <f {5 > depends, on , albeit weakly, the calcu- 

lations can, and perhaps should, be repeated a posteriori 
if more precise values of Tcj> become available. 

We now return to the task of forming a suitable sym- 
metric double difference. We shall use 



where we defer to the following section a demonstration of 
the efficacy of this definition. (Here the superscript 2s 
denotes symmetric double difference.) From Eqs. (l)-(4) / 


(6) , and (11) we see immediately that 




ft..., t)sc Tl% i lU)-%( t ...) t i'-O-'i,, 

{it) 


n ^r /if ) - Ujj ; , 
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where 


=■ fc; P U*, £*- r; (£,)). 


,|5 C p '"*1 Hj? ^ / * 0 ^ 

To discuss the implications of Eq. (12) conveniently, we 

shall introduce some approximations. First we introduce 
— > 

the vector As : 

/ b' 


(°v ~ pf ~h > 


Cn) 


which, if the time arguments coincide, represents the vector 
separation of probes ^ and with extending from ^ to p m 
Suppressing time arguments for simplicity, we can expand 

OfJ 


in terms of R, : 

ip 


— > 


whence 


S ^ ht ’ 

= ^ ^ A r fu ) 

~ £i t + A4jr) > 

where X =(x/ixi) signifies a unit vector. The neglected 
terms in Eq. (16) will in magnitude always be less than 


Ql) 




Using Eq. (16) in Eq. (12) yields 
2 0O\ /?(«-) 21 
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By use of further approximations, the expressions' for 
the unit vector differences can be made more ^rsf>icnot\s. 

Thus, 

— > 

^ tT r A . _ ir. 

0*) 


& _ ^ f /V~ rf 

Ki r ' ’ 


and 


I Zy *fi ~f]C [ = \ r t ;-ir e/ ^Z-? r ) Hrrfryfi-fi)] 




whence 




if - 1 + (^’> i ■ 

— ^ ^ * £ *cv- x t “ /V 3 ) ^ ^ z 

where in the las't line we made use of the vector identity: 

3 * (t c) = (a. c) t - (Z-Z) C . ( 2 .,) 

The desired expression for the unit vector difference is 
therefore 

where 

— 9 — > 

Ar Vj 5 rj - r- ( 2 ^ 

is the baseline vector extending from tracking station i 
to tracking station j_. 
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What is the physical interpretation of the vector 

triple product appearing in Eq. (22)? It is simply the 

7 > 

vector obtained by projecting ,{-• onto the plane normal to 

A 

. In terms of the baseline vectors, we find by sub- 
stitution of Eq. (22) into Eq. (17): 


— ■ >c*j r> , 


t!r 

(J«j 


where we also neglect the very small difference between r^ n ^ 

(n-1) a. — > ev 

and r ev * So lon 9 a s the vectors ) for the 

different pairs (ij) are not parallel, Eq. (24) will allow 
the changes in the vector separation of the probes, projected 
onto the plane normal to the earth-Venus line, to be followed 
during the period of continuous tracking of the probe pairs. 
The actual projected vector separation, as opposed to changes 
in it, cannot be determined from these data alone because the 
initial such separation — at the time simultaneous tracking 
commences is uncertain due to the fringe ambiguity. Because 
of the narrow band of the emissions from the probes, the pro- 
jected separation of the probes will only be determined to 
within the equivalent of an integral number of fringes. A 
single fringe corresponds, for the typical VLBI baselines 
under consideration, to a projected distance at Venus of about 
3 km. This ambiguity can be eliminated by simultaneous use of 
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a shorter baseline interferometer pair for which the fringe 
spacings in the two orthogonal directions, in the plane normal 
to the earth-Venus line, are larger than the corresponding 
a P r ^ or i uncertainties. This elimination is useful for 
establishing the geometry. 

In particular, our main object is to determine wind speeds 

To this end, we develop the expressions for the components of 

the velocity difference on the plane normal to r &v 

(a dot signifies differentiation with respect to time) . For 

convenience, we ignore the superscript (n) , assume we have a 

continuous determination of the projection of />,■ ; p (i) , and 

/ V 

define 

T' ? J- . „ , ^-1 ANA - A \ 'V r 


A « * (/$■*.) i * (/»•*>* 


where the unit vectors and e 2 are mutually orthogonal and 
lie in the plane normal to r ev * Since the relevant portion 
of the probes' descent through Venus' atmosphere occurs on 
a time scale short compared to a day, we ignore here the time 
dependences of b- • and * ev «* The vectors <3 W C^'b z ) can be 
defined, for example, by 


£ - ^ x 

^Li — 


f- 






A 


A 

i *f f i 


(2fcJ 


(.2- iJ 


A 

where Cj>, is parallel to the projection of and is 


♦These and other mathematical approximations would, of course, 
not be made in a realistic model to be used in an actual 
analysis cr dati. 
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normal to this projection. We also ignore the slight 

-o 

change in direction of p during final descent insofar 
as this change affects <3^ . Both components of 
are determined since we assume that the projections of 
the various baseline vectors: span the plane normal 

to ^ . 



ev 

How may we estimate wind speeds from this measured 
vector function ? First, we assume that the horiz- 
ontal velocity of each probe is equal to that of the local 
wind (see 12") . Second, we decompose the velocity 


of each probe into its vertical and horizontal parts and 


project each onto the plane normal to to determine 


their 

effects on the measured vector function. 

, Thus we 

set 



(zS) 

where 

?y = 

(ft -ft ) A > 



A" * 

x /r) i 

C 30) 

and where we again ignore the variation in j^p 

during 


. descent in its effect on the unit vector fip . in terms 
of these definitions and similar ones for the qth probe, 


we have: 
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Prom this general expression, we can- examine some special 
cases : 


(1) One of the objects being tracked has a trans- 

ponder and is not passing through the lower 
atmosphere. Since the trajectory for this 
object, say c[, can be reconstructed from the 
Doppler data, the coefficient of &p- L can be 
isolated and will yield directly the projection 
o f . 

(2) Both objects being tracked are passing through 


the lower atmosphere, but the unit vector /», 

/ t 

and the velocity component . parallel to? are 

e v 

known for each probe from transponder data (or, 
equivalently — if it were possible, from suffi- 
ciently accurate a priori knowledge of the trans- 
mitter frequency for each >probe) . In this instance, 
the differential VLBI data can be used to yield 
the time dependence of the vector difference 
between the projections of the horizontal velocity 
components of the probes onto the plane normal to 

This function can be compared with the corres- 
ponding estimate obtained from the trajectory re- 
construction. _ 


-.(3) No data other than the differential VLBI data are 
available. Here, there are a number of subcases 
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that should be mentioned. First, assume A 

^ * < 
anc * r ev are coplanar . We would then have 

/V 


~ and - i . The com- 
bined coefficient of t and in 

r 6 

Eq. (31) would therefore yield directly the 
projection of the horizontal velocity differ- 
ence. One can then apportion horizontal vel- 
ocities between the probes in a variety of 
ways consistent with the measured function 

and with "plausibility." if S», and.? 

J i* / i 3 ev 

are not coplanar, the differences of the pro- 
jected horizontal velocity components do not 
separate. Several alternatives then exist: 

(i) admit defeat; (ii) arrange to have one 

of the probes enter at the subearth point so 

joV 

that f>^ would have no component in the plane 
normal to r ev / thus allowing the coefficient 
to depend only on the projections of 
horizontal velocity components; or (iii) use 


all other available data to estimate p ’ and 
so that the observed function can 

be used to delimit the differences in the pro- 
jected horizontal velocities.- 
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( 4 ) Null results are obtained- Suppose we get a 

null value for the difference in projected hor- 
izontal velocities. What other possibilities, 
aside from the absence of winds, would be con- 

— > i 

sistent with such a result? Unless the ^ A 
are normal to r Qv — a very unlikely event-- 
the only other possibilities are either a can^- 
cellation of the projections from the two probes 

—T ^ 

or an alignment of the projections of and 

**■*& V 

p for each probe and a consequent apparent 

1 H 

The cancellation might come 


H 

absence of 

about, for example, if the winds were east-west 
at a constant speed, independent of height, and 
if the probes entered symmetrically about the 
meridian of the subearth point - 


In summary, we have developed a mathematical model to 
show that in most, circumstances the differential VLBI meas- 
urements will yield information on the wind speeds in 
Venus' lower atmosphere. But the crucial question concerns 
. whether or not such information is useful, 
is addressed in the following section. 


That question 
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IV. Error Analysis 

For a proper assessment of this application of differ- 
ential VLBI, we must investigate a large number of possible 
sources of error. It has not been possible in the limited 
time available to carry out as complete an investigation 
as is required. Thus, we shall simply list many of the 
questions which need answering, followed in turn by the 
status of our analysis of each: 

1. What is the basic resolution capability of 
differential VLBI with respect to the entry 
probes? 

2. What are the limitations imposed by: 

i. lack of clock synchronization be- 
tween the various receiving sites 
on earth; 

ii. the receiver systems? 

iii. the atmosphere of Venus; 

iv. the atmosphere of the earth; 

v. the ionospheres of Venus and the 

earth, and the interplanetary medium? 

vi. instability of the transmitter freq- 
uency? 

vii. uncertainty in^ the geometry of entry 

for the probes or in the trajectory of 
the biAJ* if the latter is used as a 


reference? 
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3. What are the optimum configurations of the 
probes with respect to positions and times 
of entry? 

1. To determine the basic resolution capability, we 
note that the weakest signals will be from the miniprobes 
which will transmit about 1 w of effective radiated power 
when near the surface of Venus. If the bandwidth of this 
signal is no more than 50 Hz (probably a gross upper bound), 
then the flux at the earth will be no less than about 30 FU 
(1 FU = 10 w/m -Hz). For the antenna systems that might 
be used in the experiment — Golds tone , Madrid, 

Arecibo, Haystack, and Johannesburg--the fringe phase un- 
certainty, due solely to system noise, would be under 1 ° 
after only one or two seconds of integration (see, also, 
2.ii below) . Such a fringe phase error corresponds to a 
displacement uncertainty at Venus of about 6 m for this 
S-band signal with a projected baseline of 4000 km and 
an earth-Venus separation of 0.5 a.u. Thus average pro- 
jected velocity differences could be measured over a time 
•interval t with an uncertainty of only about 10 m/sec, 
where t is in seconds, if the system noise were the only 
source of uncertainty. We may compare this resolution 
with the time required for the entry probes to acquire the 
horizontal speed of the wind. As a crude model, con- 

sider the probe to be spherical of radius R and average 
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density y9 . Then, if we neglect the density of the at- 
mosphere relative to that of the probe, the probe's hor- 
izontal acceleration a R will obey: 


*h * * -/*U) , 

(3iJ 

whence its velocity v H . will be given by 


v l e / 


where the time constant X” 1 is 


K -'= , 

Cs 1 *) 

with \ being the viscosity. Since the 
velocity 'ly is given approximately by 

terminal vertical 

av - > 

(3$) 

we have X ^ 5 sec for ^\T y - 50m/sec, 

etc* The required 


VLBI integration time thus appears well matched to the time 

scales in which the probes reflect the local wind speeds in 
the lower atmosphere. 

2.i. The lack of precise clock synchronization between 
the various receiving sites should introduce no detectable 
error if the data are properly taken. With the signals from 
each object tracked being sampled simultaneously, the clock 
error cancels completely upon differencing. In effect, the 
strong signal, say from the bus or main probe, acts as the 
clock for the weaker signals from the miniprobes. A large 
epoch offset of the station clock from one site relative to 
that from another only increases the set of trial times that 

* Here 2 is the acceleration of gravity c m the surface of 
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need be introduced in the usual cross-correlation procedure 
used to search for fringes. However, even this minor problem 
disappears if advantage is taken of the presence of the 
signal from each source as we explain below. 

2.ii. The receiver system, if properly configured, 
also need introduce 1 no* detectable errors. We require here 
that the different local-oscillator signals, used at a given 
site to heterodyne the radio-frequency signals from the 
various objects being tracked, all be derived from the same 
frequency standard and, insofar as possible, from the same 
L.O. chain elements. The purpose of these strictures is 
to insure that almost all of the phase noise of the hetero- 
dyne signals are common to the receiver chains for all tracked 
objects. The common phase noise thus introduced will then 
cancel upon differencing. The residual (non common) phase 
noise can probably be reduced without much difficulty to 
the order of 1°. 

In connection with the receiver system, we also note 
that the presence of a carrier signal— lacking in the usual 
* celestial sources involved in VLBI experiments—allows the 
tape recording of the heterodyned signals and subsequent 
cross correlation to be eliminated. If the carrier signal 
from each probe is sufficiently stable, it can be tracked 
with a suitable phase-locked loop (of third order, if necessary. 
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to follow drifts in transmitter frequency) and only the 
usual counted-Doppler values need be recorded. These 
samples, can be incorporated directly into the double-dif- 
ference observable defined in Section III. We must still 
insure that the samples for each tracked object are obtained 
simultaneously, or very nearly so, to insure that the clock 
synchronization errors cancel. 

2.iii. The atmosphere of Venus can be expected to in- 
troduce sizable phase variations in the signals received 
at a given site from a given tracked object. The one-way 
electrical path length of the Venus atmosphere is about 
300 m in the zenith direction. However, the phase var- 
iations introduced will be virtually identical in their 
effects on the signals received at each of the earth-based 
tracking stations. The geometric beams from a given entry 
probe to each of the tracking stations are separated by 
about 1 m at an altitude of 20 km. The Fresnel zone at 
that altitude for these S-band signals measures about 70 m 
across. Since the overlap is almost complete, this error 
source will largely vanish in the symmetric differencing 
process. The residual phase noise will be due to the small 
crescent-shaped non-overlap regions, separated by about 70 m 
at a 20 km altitude and by less at lower altitudes. This 
noise will depend on the spatial spectrum of the atmospheric 
inhomogeneities and on the wind speeds. No attempt has yet 
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been made to estimate this contribution quantitatively. 

2. iv. The earth's atmosphere introduces far less 
severe phase fluctuations since it has a zenith electrical 
path length of less than 3m — two orders of magnitude smaller 
than for Venus. The almost, complete overlap of the beams 
entering a given antenna aperture from the various entry 
probes insures a high-order of cancellation of the atmos- 
pheric effects in the differencing procedure. No quanti- 
tative estimate has yet been made of the residual noise, 
although relevant noise statistics are available. In 
summary, the Venus atmospheric effects tend to cancel 
because of the differencing of the signals received at 
the different receiving sites whereas the corresponding 
effects of the earth's atmosphere tend to cancel because 
of the differencing of the signals from the various probes. 
The beauty of the symmetric double-difference technique 
is thus apparent. 

2.v. The effect of the charged particles along the 
propagation paths — equivalent to a change in electrical 
• p&th length of less than 15 m— — will also tend to cancel 
in the double difference. But here there are several im- 
portant differences from the atmosphere case: (a) The 

path separations midway between earth and Venus are about 
2000 km, compared to a Fresnel zone size of about 100 km; 

(b) The ionospheres have peak densities at altitudes of 
hundreds of kilometers; and (c) The lack of exact equality 
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among the transmitter frequencies will prevent complete 
cancellation of plasma effects due to dispersion. To 
insure that the frequency differences A - ? cause corres- 
ponding variations in phase path of no more than 1° f it 
is necessary that Af & if the charged-particle 

contributions can be modelled from other data to within 
4 0%, then this restriction can be relaxed to Af 6 S~O.UH^ 
The lack of cancellation due to non-overlap of the various 
paths again will depend on the spatial and temporal spectra 
of the inhomogeneities. And again no attempt has yet been 
made to estimate quantitatively this source of residual 
phase error. 

It may, in fact, be possible to solve for the above 
propagation medium effects if all probes can be tracked 
simultaneously from more than three earth-based antenna 
sites. The multiplicity of paths provides redundancy which 
may be used in a suitable filtering scheme to eliminate all 
medium effects. (We assume that the differences in trans- 
mitter frequencies introduce negligible dispersion.) The 
• analysis of this multi-probed many-sited situation has been 
started, but not completed. 

2.vi. Variations in the frequencies of the transmitters 
make difficult their calibration by means of the one-way 
Doppler values. Errors in this calibration will tend to 
introduce asymmetries into the double-difference observable 
with the consequence, for example, that the propagation 
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medium effects will not cancel as completely. An a priori 
knowledge of these frequencies, accurate to 1 part in 10 8 , 
would certainly be sufficient to eliminate this source of 

error. Whether or not a knowledge of them to 1 part in 

__6 

10 --a more realistic figure — will, be sufficient is un- 
certain. A detailed analysis of this aspect is in progress 
but has not been completed. With turn-around transponders , 
the problem all but disappears. 

2.vii. Uncertainties in. the geometrical configuration 

~f> 

of the entry probe vectors P ? (see Section III) will 
affect the interpretation of the data in terms of wind 
speeds. Similarly, uncertainties in the velocity vector 
of the bus, if it is used as a reference, will introduce 
interpretation difficulties. Although no quantitative 
estimates have been made, the trajectory reconstruction 
for the bus should be sufficiently accurate with Doppler 
errors at the 1 mm/sec level or below. Hopefully, the 
medium effects and the unknown harmonics of Venus * grav- 
itational field won't vitiate this conclusion. With the 
• bus providing a reference, the small uncertainties in the 
entry probe geometry will be of little consequence? such 
will not be the case if only the several miniprobes are 
tracked simultaneously. But then other problems loom larger, 
as mentioned in Section III. 
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3- It would be best to have the bus tracked simul- 
taneously with each entry probe (no entry occultations 
allowed during this period!) and to have the trajectory 
of the bus passing far enough from Venus to minimize the 
effects of the higher harmonics of the gravity field. 

against this requirement will be the reduction 
in cancellation of propagation medium effects that accomp- 
anies an increase in angular separation of the targets, 
if the bus is unavailable, then wind speeds seem to be 
easiest to isolate (see Section III) if one of the simul- 
taneously tracked entry probes is directed towards the 
subearth point on Venus. Again, the quantitative advantages 
have not been analyzed. 
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V. Conclusions 

We conclude that wind speeds in the lower atmosphere 
of Venus can be detected via differential VLBI observations 
of the entry probes;. The uncertainty in the wind speed 
determination can probably be kept below 100 t~^m/sec, 
where t > 1 is the time resolution in seconds, provided that: 

(1) the residual effect of the propagation medium on 
the symmetric double-difference observable can be 
be kept below about 10° of phase at S-band? and 

(2) the transponded signals from the bus, on an exo- 
atmospheric trajectory, are available as a reference. 

Condition (1) , which is crucial, is unfortunately not but- 
tressed by a prima facie case. The loopholes left, discussed 
in Section XV, are related to the effects of the non-overlap 
regions of the propagation paths and the instabilities of the 
transmitters. (We assume, in addition, that the differences 
in transmitter frequencies are no more than about 50 kHz .) 

If the signals from the bus were not available as a reference, 
the main probe with its transponder could serve the same 
* function with a loss in accuracy that would probably not be 
too severe but that hasn't been estimated quantitatively. 

If only the miniprobes — without transponders — can be tracked 
simultaneously, the situation looks grim because of the 
difficulty in separating the contributions of the vertical 
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and horizontal velocity components to the observed pro- 
jections on the plane normal to the earth-Venus line. 

Unless the vertical velocity components can be modelled 
accurately, the only apparent solution in this circumstance 
is to have one of the miniprobes directed towards the 
subearth point. 

Of course, if a transponder were available on each 
probe, all of these problems would fade away. The recon- 
struction of the horizontal probe velocity from the pro- 
jection of the total velocity vector along the earth-Venus 
line should be reasonably accurate even after allowance for 
uncertainties in the entrance geometry, the terminal vertical 
velocity, and the (two-way) effects of Venus' atmosphere on 
the observed Doppler shift. 

The main conclusion to be drawn on the potential of 
differential VLSI, per se, for the determination of wind 
speeds is that, despite this mass (mess?) of verbiage, much 
work remains to be done to assure a proper assessment. 
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Astronomical Applications of Differential Interferometry 

Abstract. Intercomparison of radio signals received simultaneously at several 
sites from several sources with small mutual angular separation provides a power- 
ful astrometric tool. Applications include tracking the Lunar Rover relative to 
the Lunar Module, determining the moon’s libration, measuring winds in Venus’s 
lower atmosphere, mapping Mars radiometrically, and locating the planetary 
system in an inertial frame. 


In most applications of very-long- 
baseline interferometry (VLBI) the 
most serious limitations on the accuracy 
of the results are imposed by unknown, 
variable phase errors introduced by both 
the neutral atmosphere and the iono- 
sphere above the receiving sites, and by 
fluctuations in the rates of the oscillators 
that provide phase references at the 
separate sites. These limitations may be 
largely removed in differential measure- 
ments, in which signals received simul- 
taneously from different radio sources 
located close together in the sky are 
compared. If atmospheric and indepen- 
dent oscillator phase shifts affect obser- 
vations of each source equally, their 
effects will cancel when differences be- 
tween observations are examined. In this 
report we discuss several scientific ap- 
plications of differential interferometry 
(i), as well as the actual tracking of 
the Lunar Rover performed during the 
Apollo 16 mission. 

Because differential interferometry 
involves taking differences not only be- 
tween receiving points but also between 
transmitting points, it follows that any 
potential source of error will cancel if 
it is common either to all receivers or 
to all transmitters. This simple principle 
will be shown to have important con- 
sequences for astronomical measure- 
ments. One such consequence relates to 
observations of artificial transmitters for 


which the carrier frequency may be un- 
certain and variable. Noninterferomet- 
ric one-way Doppler tracking of such 
objects is ordinarily of little use be- 
cause changes in the received frequency 
due to the Doppler shift cannot be 
distinguished from changes in the fre- 
quency of the transmitter itself. In inter- 
ferometry, however, transmitter (and 
any other) frequency changes that ap- 
pear equally at all receivers have no 
direct effect on the ability of the inter- 
ferometer to determine relative angular 
positions. In faot, artificial radio sources 
make particularly convenient objects 
for interferometry because conventional 
Doppler counting techniques can be 
used to keep track of the phase of 
the carrier signal received at each site. 
Wide-bandwidth group-delay interfer- 
ometry also may be done efficiently 
with artificial sources if the carrier wave 
is suitably modulated, for example, with 
a pseudorandom wave form of the kind 
often employed for two-way radar rang- 
ing (2). As in the case of one-way 
Doppler tracking, one-way radar rang- 
ing is ordinarily useless if either or 
both of the transmitter and receiver 
time bases are unstable. But for either 
phase-delay or group-delay observables, 
the effects of transmitter instability can- 
cel when the difference is taken be- 
tween receiving sites (thus forming an 
interferometric observable), and the ef- 


fects of receiver instability cancel when 
the difference is taken between a pair of 
transmitters (forming a differential in- 
terferometric observable). 

We shall now discuss some of the 
potential scientific applications of dif- 
ferential interferometry. First, however, 
we describe one technical application 
already successfully carried out: Earth- 
based tracking of the Apollo 16 Lunar 
Rover relative to the Lunar Module. 
Three tracking stations (3) were em- 
ployed so that two independent base- 
lines were formed. Thus, two compo- 
nents of the motion of the Rover rela- 
tive to the Module were determined 
from the changes in phase of the two 
differential interferometric observables. 
From the initial separation of the Rover 
and the Module and the constraint that 
the Rover remained on the lunar sur- 
face, it was possible to determine its 
entire path (Fig. 1). After a traverse of 
over 4 km, the final position computed 
from these data differed from the actual 
position by about 30 m, or about 0.01 5 
arc second at the distance of the moon 
(4). The main source of error was rela- 
tive phase drift between the two re- 
ceivers (one each for the Rover and 
the Module) used at each site in this 
trial experiment. In an operational sys- 
tem this error would be eliminated by 
using a single receiver for both signals. 
The basic technique appears capable of 
reducing tracking errors to the meter 
level, a limit imposed by unmodeled 
lunar topography. 

A related scientific application in- 
volves the accurate determination of the 
moon’s libration by monitoring simul- 
taneously from several tracking stations 
the ALSEP <5) telemetry transmitters 
located at three well-separated sites, 
such as those of Apollo 14, Apollo 15, 
and Apollo 16. Here, because the 
ALSEP’s are fixed on the lunar surface, 



Lunar 

Module 



Fig. 1. The path of the Apollo 16 Lunar Rover is shown as 
determined by Earth-based differential interferometric tracking 
on 21 April 1972. Individual dots mark the positions obtained 
at 20-second intervals, beginning at 20:52:40 U.T. from point 
A. Craters given names by the astronauts are included for refer- 
ence, although their locations are known only approximately. 
The Rover was stopped at point B for 6 m 20% at C for l h 9 ,n 40", 
and at D for 27 m 0'*; several brief stops were made at E. At 
23:03:40 our tracking indicated that the Rover had stopped 
finally at F, 30 m east of the Lunar Module; the Rover had 
actually parked at the Module. Some of this error may reflect 
a corresponding error in the assumed starting position, A. How- 
ever, tracking data obtained while the Rover was known to 
be stopped occasionally showed systematic drifts as large as 
2 or 3 cm/sec. (see text). Random noise was less than 1 m. 
At all times during the traverse, position readings "from the 
navigation system on board the Rover agreed within 100 m 
(approximately the limit of precision of the onboard system) 
with these differential interferometric tracking results. 



their relative positions must be deter- 
mined by monitoring changes in the 
phases of the differential interferometric 
observables over a sizable fraction of a 
day. Lunar libratiop causes these ap- 
parent positions to vary over the course 
of a month, and longer. The differential 
nature of the observable sharply reduces 
the effects of errors in the lunar ephem- 
eris, tracking station coordinates, and 
so forth, and should yield at least an 
order of magnitude improvement in our 
knowledge of the libra tion [present un- 
certainty about 10 seconds of selenocen- 
tric arc (6)]. 

Differential interferometric tracking 
of planetary probes, landers, and orbit- 
ers will yield results in many cases more 
accurate than can be obtained from 
tracking satellites of the earth. This 
seemingly paradoxical conclusion fol- 
lows because the usual limitation is set 
not by signal strength but by systematic 
effects of the atmosphere and iono- 
sphere, and sometimes by receiving- 
system phase instabilities. These effects 
cancel in the differential interferometric 
observable. A planetary application of 
differential interferometry which is 
analogous to, but more complicated 
than, the Lunar Rover-Lunar Module 
situation involves tracking a number of 
small probes descending simultaneously 
into Venus's atmosphere (7). By differen- 
tial tracking of the free-falling probes 
relative to a parent spacecraft it should 
be possible to detect horizontal winds at 
the level of a few meters per second. 
Differential interferometry could also 
aid in the interpretation of occultation 
data (5). Additional applications of dif- 
ferential interferometry to both orbiters 
and landers are too numerous to be 
elaborated here; for example, improved 
estimates of the planet's gravity model 
parameters, rotation vector, and landing 
site or orbit parameters may be ex- 
pected [see, for example, (5)]. Whenever 
the planet passes close to the direction 
of an extragalactic radio source, dif- 
ferential interferometry may be used to 
determine the earth-planet direction 
relative to that of the source to about 
0.001 arc second. Such measurements 
could be used to determine precisely the 
orientation of planetary orbits with re- 
spect to an inertial frame and, for ex- 
ample, to monitor the perihelion preces- 
sions to further test general relativity. 
For ground-based radiometric mapping 
of the terrestrial planets, differential in- 
terferometry can overcome the effects 
of instrumental, atmospheric, and iono- 
spheric phase drifts, which limit the ap- 
plication of aperture-synthesis tech- 


niques. A phase reference could be pro- 
vided, for example, by the specular echo 
of a radar signal sent at the appro- 
priate. frequency (9). Such mapping ap- 
pears especially important for Mars 
where the distribution of small amounts 
of surface water (or ice) might be dis- 
cernible from millimeter-wavelength ob- 
servations (10). 

In summary, the technique of differ- 
ential interferometry seems capable of 
solving a wide range of astronomical 
problems. 

C. C. COUNSELMAN, III 
H. F. Hinteregger 
I. I. Shapiro 

Department of Earth and Planetary 
Sciences, Massachusetts Institute of 
Technology, Cambridge 02139 
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Venus; Radar Determination of Gravity Potential 


is given approximately by the semi- 
empirical formula (6) 


Abstract. We describe a method for the determination of the gravity potential 
of Venus from multiple-frequency radar measurements. The method is based on 
the strong frequency dependence of the absorption of radio waves in Venus’ 
atmosphere. Comparison of the differing radar reflection intensities at several 
frequencies yields the height of the surface relative to a reference pressure 
contour; combination with measurements of round-trip echo delays allows the 
pressure, and hence the gravity potential contour, to be mapped relative to the 
mean planet radius. Since calibration data from other frequencies are unavailable, 
the absorption-sensitive Haystack Observatory data have been analyzed under the 
assumption of uniform surface reflectivity to yield a gravity equipotential contour 
for the equatorial region and a tentative upper bound of 6 X 10~'> on the frac- 
tional difference of Venus’ principal equatorial moments of inertia. The minima 
in the equipotential contours appear to be associated with topographic minima. 


Present knowledge of the surface of 
Venus rests largely on the results of 
radar observations. Perhaps the most 
striking fact to emerge has been the 
retrograde direction of Venus’ spin and 
its apparent resonance with the relative 
orbital motions of the earth and Venus 
(/). The earth could have captured 

Venus' spin in this resonance only 
through the action of a gravitational 
torque on a substantial axial asymmetry 
in Venus’ mass distribution (2). Hereto- 
fore, no measurement of this asymmetry 
has been possible. The main purpose 
of this report is to demonstrate that 
future radar observations can be used 

to determine equipotential contours of 
Venus' gravity field and, hence, to 
estimate the axial asymmetry of its mass 
distribution. A preliminary contour for 
the equatorial region and a concomitant 
bound on the axial mass asymmetry — 
based on past radar observations not 
made explicitly for this purpose — is also 
included. 

How can radar data be sensitive to 
the gravity field of Venus? A direct 
sensitivity seems almost unthinkable. 
But an indirect intermediary exists, 
namely, the thick, carbon-dioxide-domi- 
nated atmosphere of Venus. Because 
this atmosphere absorbs X-band (ap- 
proximately 8000 Mhz) radio radiation 
strongly and, for example, S-band (ap- 
proximately 2000 Mhz) radiation hardly 
at all, we can infer surface heights 
relative to a particular pressure contour 
from a comparison of radar cross sec- 
tions measured at the two frequencies, 
since the intrinsic reflectivity of the 
surface itself should not, in general, 
vary sharply with frequency (3). The 
use of a third frequency would allow 
a more precise separation of atmo- 
spheric from surface reflectivity effects 
on cross section. The pressure contours 
can then be related to the mean planet 


radius with the aid of measurements of 
round-trip radar echo time-delays, 
which allow the absolute surface heights 
to be determined (4). Gravity equipo- 
tential contours will coincide with pres- 
sure contours under conditions o‘f hy- 
drostatic equilibrium in the atmosphere 
(5). From such a contour, the gravita- 
tional torque exerted by the earth can 
be estimated.. 

Now we develop this basic idea 
quantitatively. Since at present its ap- 
plication is restricted to the equatorial 
regions traversed by the subradar point 
of earth-based observations, we confine 
our analysis to that situation. The pos- 
sibilities for extension to high latitudes 
and for the use of radarj on Venus 
orbiters are discussed briefly in the last 
part of the report. 

The radar cross section <r(A,0) per 
unit surface area at the subradar point 
can be written as 

a(\,<p) = OY,(X,0)exp(— 2r(\,4>)] (1 ) 

where 

t(X,0) := J «(.h,\)dh (2) 

H<P) 

with A being the wavelength of the 
radar signals, 0 the longitude of the 
subradar point (we suppress 6, the lati- 
tude dependence), the intrinsic cross 

section per unit area of the observed 
region of the surface, r the opacity 
(optical depth) of the atmosphere, k the 
absorption coefficient for radio waves, 
h the height of the reflecting region 
relative to a reference pressure contour, 
and h [n&x the altitude above which ab- 
sorption can be neglected. We define 
the reference contour in terms of a 
reference longitude 0 O such that 
A(0 O ) 23 0. The factor of 2 multiplying 
r in Eq. 1 accounts for the two-way 
passage of the radio waves through 
Venus’ atmosphere. The coefficient, k. 
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where P is the pressure in atmospheres, 
T the temperature in degrees Kelvin, 
and A the wavelength in centimeters. 
Since the absorption is important only 
in the lower atmosphere, we may use 
the approximate temperature-pressure 
relation (7) 


PCA)^r r(O) -Lh y 

P( 0) L T(0) J 1 J 

in the evaluation of t. Here h = 0 refers 
to the reference pressure contour; L ~ 
9°K per kilometer is the lapse rate, 
assumed constant' and k = gg/RL ^ 5.3 
is the polytropic index, with p. the mean 
molecular weight of the atmosphere, g 
the acceleration of gravity, and R the 
gas constant. From Eqs. 1 through 4 
we obtain 


g(A .0) tro(X,0o) ^ 
ff(X,0 o ) <r o(X,0) 

Mr<0) 

2 )L 


exp 
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« exp[2*(0,X)/i] (5) 

where the last relation is valid only for 
h « [T{0)/L\ * 90 km. Since h does 
not seem to vary by more than about 
± 4 km ( 4, 8), the expansion used in 
the last part of Eq. 5 will be in error 
by less than 15 percent. Hence, for 
expository purposes, we confine discus- 
sion to the simplified form. Solving for 
h then yields 


A(0)~ 


t . r g(x,-») <7 „(x,«ft„) ~i ... 

2«(O,X) Jn U(X,0 o ) *,(M)J ^ 


Assuming that the multiple-frequency 
observations provide the ratios a! <j 0 , 
how may we use the resultant values of 
/j(0) to determine the height variation 
of the reference contour with respect 
to the mean surface radius, p? From 
the value of the echo time-delay, mea- 
sured simultaneously with the cross sec- 
tion, we can infer the height, h f , of the 
reflecting region above the mean radius 
(4). Thus, the height, h", of the 
reference pressure contour above the 
mean radius is given by 


h’’ = h'-h (7) 

The function defines a gravity 

equipotential contour over the equa- 
torial region under conditions of atmo- 
spheric hydrostatic equilibrium (5). If 
the contribution of the centrifugal ac- 
celeration to the gravity field were 
neglected (9) and if h n were known 
over the entire planet, then the un- 
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known coefficients in the expression for 
gravitational potential energy, U , could 
be obtained from inversion of 

CMg 

u[p + = 7-rx x 

p + n 

(' + E E(7f7r)' p ’< sta8 ) x 

\ # = l n =o' fT ' 

= C/». (8) 

where G is the gravitational constant, 
M s the mass of Venus, ^(sin^) the 
associated Legendre function of degree 
n and order m, C„ m and S„ M the sought- 
for coefficients, and C/ 0 the value of the 
potential on the reference contour. For 
data confined to the equatorial regions, 
as here, Eq. 8 can be recast as 

( C,„cosm<fi -f S„,sinm<t>) ~ ^ ^ (9) 

m=-n p 

where 


fC,,mCos/?i0 -f S nm smffifp] 



arid where we set U 0 =* — GM^/p (10) 
and dropped terms of second order and 
higher in h The coefficients C,„ and 
S,„ are given by 

2 7T 

r^l=-L f h"{<P)[ cosm f]d<p; 

LS» J *p Lsinm^J 

m = 1 , 2 ,... ( 11 ) 

with C 0 ^ C 00 = 1 and = 0. No use- 
ful information on the coefficients of 
the zonal harmonics (C n0 ; « ^ 2) is con- 
tained in the equatorial portion of the 
equipotential contour, since the zonals 
have no longitude dependence and their 
bulk equatorial effect is not easily 
separable from that of 

How may We use these results to esti- 
mate the gravitational torque exerted 
on Venus by the earth? For study of the 
putative spin-orbit resonance (2), the 
relevant torque is proportional to (B — 
A)/C where A < B < C are the princi- 
pal moments of inertia of Venus, with C 
assumed to be the moment about the 
spin axis. Unless either (i) Venus is now 
a very elastic body (high “£>’’) with 
respect to the diurnal stresses of 100- 
day periodicity, or (ii) there exists a 
very delicate balance between the tor- 
ques exerted by the sun on Venus’ tidal 
bulge and on a possible atmospheric 
bulge (//), it appears that control of 
Venus’ spin by the earth requires 
(B — A)/C > 10— 1 (2). If the tesseral 


harmonics for Venus fall off with degree , 
as do those for the moon and Mars, 
then the approximation 

i(c*y+ <&)*]* ^ 

3[(C m ) 2 + = 0.3^^ (12) 

should yield a reasonable estimate for .. 
(B-A)/C. 

Unfortunately, data necessary to de- 
termine an accurate equipotential con- 
tour in the equatorial region of Venus 
do not now exist. The lack of both 
accurate values for *:(h,\) and properly 
calibrated radar cross-section data are 
the major limitations; in particular, there 
have been no coordinated observations 
of cross section at more than one radar 
frequency. If the variations with longi- 
tude of the intrinsic surface reflectivity 
are small and if the effects of the differ- 
ences between total radar cross sections 
and cross sections per unit surface area 
at the subradar point are also small, we 
may use the limited data from the 
Haystack Observatory on the total cross 
sections (8) at X = 3.8 cm, coupled with 
the surface-height variations recently de- 
termined (4), to obtain an approximate 
equipotential contour. The result is pre- 
sented in Fig. 1 for k = 0.07 km -1 . 
This choice for k is based on Eq. 3 and 
the “nominal” values P(0) = 100 atm 
and T(0) = 750°K (12). A comparison 
with the surface-height variation (see 
Fig. 1) seems to indicate that minima 
in the equipotential contour are as- 
sociated with topographic minima. But 
one must remember that the uncertain- 
ties are large; it is even difficult to place 
reliable bounds on the accuracy of the 
equipotential contour, in view of the 
lack of accurate and suitable data. 

We can assess analytically the rela- 
tive sensitivity of h " to the various 
relevant factors. From Eqs. 1 through 
7 we find, under conditions validating 
the last part of Eq. 5, 


W = 8h' + 


1 5g 
2«(0,X j a 


Sk( 0.X) 
40, X) 


(13) 


The last term, through Eq. 3, can be 
separated into components due to errors 
in 7X0), P(0), and numerical factors. 
On the basis of experimental and 
theoretical evidence relating to these 
components, we estimate that \8k/k\ 
<0.6. This contribution to h” is pro- 
portional to h and may therefore be 
as much as several kilometers. For our 
nominal value of «(0,A = 3.8 cm) — 0.07 
km" 1 , fractional errors in cross section 
of 10 percent will contribute errors of 
under I km to h". The contribution of 


8 h f (4) should be nowhere greater than 
about 0.5 km. 

Unless the intrinsic surface reflec- 
tivity is well correlated with the surface 
altitude — unfortunately a not. unlikely 
possibility — the solutions for C 2 and S 2 , 
and hence the estimate for ( B — A)/C, 
may be relatively immune to the effects 
of variations in intrinsic reflectivity. Un- 
der this assumption, we find from Fig. I 
and Eqs. 11 through 13 that 

=* (3 ± 3) x 10- 4 (14) 

where the error reflects our estimates 
of the uncertainties in k, c(0), and 
h'(0). No allowance was made for the 
contributions of the higher-degree terms 
to C 3 and S. 2 (13). We also find that the 
axis of minimum moment of inertia 
passes through longitude 30 ± 60 de- 
grees [International Astronomical Union 
(IAU) coordinate system], which is too 
uncertain to allow meaningful deduc- 
tions about torque balance. . 

The prospects tor improvement in 
this very crude estimate of the gravity 
equipotential are good. Results from 
Venera 8, for example, should tighten 
the bounds on atmospheric composition, 
7X0), and P(0), and hence on k (after 
laboratory confirmation or correction of 
Eq. 3). The radar systems at Arecibo, 
Puerto Rico (X = 70 cm), Goldstone, 
California (X = 12.5 cm), and at Hay- 
stack (X = 3,8 cm) could be used 
to reduce the uncertainties in surface- 
height variations to the 150-m level or 
slightly below. With careful calibration 
of the. radars, cross section measure- 
ments — for the same surface regions to 
which the height measurements apply — 
should have relative errors of no more 
than about 2 percent (8). Haystack’s 
contribution is essential here, because 
there is no appreciable atmospheric 
absorption at the wavelengths used at 
the other observatories. Thus, earth- 
based measurements could yield gravity 
equipotential contours in the equatorial 
regions of Venus with a lateral surface 

resolution of about 100 km, correspond- 
ing to information on spherical harmon- 
ics up to the 360th degree, and an 
altitude resolution of about 200 m. The 
decrease in the uncertainty of the esti- 
mate of (B — A)/ C should be at least 
fourfold if the higher harmonics do 
not contribute too much to C 2 and S. 2 . 

Can the determination of Venus’ 
gravity equipotential contour be ex- 
tended beyond the equatorial regions 
traversed by the subearth point? Two 
approaches are possible, (i) With more 

y 

/ 

/ 

/ 



large (see text). Latitude variations have been suppressed, witn 
a narrow (± 10°) band centered on the 


Fig. 1. Comparison 
of the surface 
heights on Venus 
■with a gravity equi- 
potential contour. 

The surface heights 
are based on round- 
trip radar echo de- 
lays from published 

( 4 ) and recent, un- 
published data. The 
gravity potential 
contour is derived 
from data on radar 
cross sections ob- 
tained at the Hay- 
stack Observatory 

(5) ; the uncertain- 
ties in both the scale 
and variations of 
this latter curve are 
both curves referring to averages along 
equator of Venus. 


powerful earth-based radar systems, 
such as the proposed improved Arecibo 
facility, it will be possible to determine 
surface heights and reflectivities over 
most of the planet with high resolution 
by use of the new technique of delay- 
Doppler interferometry (14) at the 
12.5-cm wavelength at which this radar 
would operate. If a similar capability 
existed for shorter radar wavelengths, 
such as with a Haystack-Goldstone 
bistatic cqnfiguration, then the atmo- 
spheric absorption could be determined 
as well. Of course, in the analysis of 
these data — for which the incident and 
reflected waves would not be in the 
zenith direction on Venus — atmospher- 
ic refraction effects must be considered 
(7) as well as possible variations of the 
intrinsic angular scattering law with 
frequency (3), (ii) A spacecraft placed 
in a polar, or near-polar, orbit about 
Venus and equipped with a suitable 
dual-frequency radar, could determine 
surface heights, reflectivities, and the 
corresponding atmospheric absorptions 
over virtually the entire planet. Re- 
peated polar passages would offer the 
possibility for continual calibration. The 
resultant high-resolution equipotential 


contours would not only yield the 
gravitational torque exerted by the earth 
but would have other scientific applica- 
tions as well: these contours bear, for 
example, on questions of the origin 
and evolution of Venus, its deep inte- 
rior, the extent of isostatic compensa- 
tion near its surface, and the processes 
of surface erosion (15). 
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